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Abstract 

Triple differential cross sections of pions in heavy ion collisions at 1 GeV/nucl. 
are studied with the IQMD model. After discussing general properties of 
A resonance and pion production we focus on azimuthal correlations: At 
projectile- and target-rapidities we observe an anticorrelation in the in-plane 
transverse momentum between pions and protons. At cm. -rapidity, however, 
we find that high pt pions are being preferentially emitted perpendicular to 
the event-plane. We investigate the causes of those correlations and their sen- 
sitivity on the density and momentum dependence of the real and imaginary 
part of the nucleon and pion optical potential. 
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I. INTRODUCTION 



One of the main goals of the study of relativistic heavy ion collisions is the determination 
of the density and momentum dependence of the real and imaginary part of the nucleon 
(and other hadron) optical potential (often also termed nuclear equation of state) 
Its importance stretches well beyond nuclear physics and is of great importance for the 
formation of nuclear matter after the big bang, the behaviour of supernovae and neutron 
stars. It also is important for the quest for the quark-gluon plasma in heavy ion collisions. 

An increasing number of observables which are accessible through heavy ion collisions 
has been found to be sensitive to the equation of state: Among the most prominent ones 
are collective flow effects such as the bounce-off of cold spectator matter in the reaction 
plane || and the squeeze-out of hot and compressed participant matter perpendicular to 
the reaction plane || as well as particle production ||10Hl"2]1 . The pion multiplicity was one of 
the first observables suggested to be sensitive to the nuclear equation of state fl0|-|l2"||. This 
was motivating a strong experimental effort (4n analysis of streamer chamber events at the 



BEVAL AC) [|i3| -|l5ll . However, the sensitivity of pion yields and spectra [ 16 1 on the equation 



of state is not very high |17|Jlq| and therefore the attention shifted towards subthreshold 
production of mesons (e.g. kaons and 77-mesons) ||l"9|-23[ . 

New experimental 4ir setups at two of the major heavy ion reasearch facilities, GSI 
(FOPI, KaoS, TAPS) and LBL (TPC), enable the investigation of the emission pattern and 
correlations of primary and secondary particles in a far more detailed manner than ever 
before. It is now for the first time possible to thoroughly investigate correlation phenomena 
such as in-plane bounce-off P^-p7| and out-of-plane squeeze-out p8|- p0|1 of pions. The 
detailed investigation of these effects, including their possible origin and their impact pa- 
rameter and pt dependencies as well as their sensitivity to the nuclear equation of state, are 
the subject of this publication. 



II. THE IQMD-MODEL 

The first widely used microscopic models for the description of relativistic heavy ion col- 
lisions were based on the Vlasov-Uehling-Uhlenbeck (VUU) theory |T8|J31| , |32[] , which explic- 



itly treats nonequilibrium and (stochastic) quantum effects in the framework of one-particle 
quantities, as well as the nuclear potential (nuclear equation of state). The dynamical basis 
of the VUU-model is the following transport equation: 

Bf 4 r da 
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x[fh{i-f x ){i-r 2 )-f[m-f){i-f2)] 

x 5 3 (p + p 2 ~p' 1 -p' 2 ). 

f is the single-particle distribution function. The l.h.s. contains the potential U. Usually U 
is parameterized using the Skyrme ansatz. This gives the opportunity to study the effects 
of the nuclear equation of state via different parameter sets. 

The r.h.s. contains the cross section a and the Nordheim-Uehling-Uhlenbeck modifica- 
tions incorporating the Pauli-blocking factors [Q. Models based on the same theory, but 
differing in numerical implementation, are the Boltzmann-Uehling-Uhlenbeck (BUU) f32j 
and the Landau- Vlasov |34| model. All 3D-numerical implementations of the VUU theory 
are solved with the test particle method. The number of test particles used to represent a 
nucleon varies with the numerical implementation. The test particle method solvs Hamil- 
ton's equation of motion for each test particle. These transport models have been successful 
in studying various aspects of relativistic heavy ion collisions, such as single particle spectra, 
collective effects (stopping, bounce-off, squeeze-out) and meson production. 

However, certain fluctuations and correlations, such as the formation of fragments in 
relativistic heavy ion collisions, cannot be studied with a transport model based on a single- 
particle distribution function. This was one of the motivations for the developement of the 
Quantum Molecular Dynamics model (QMD) pl|j35| -p7|. In the QMD model the baryons 



are represented by Gaussian shaped density distributions 

They are initialized in a sphere of a radius R = l.lAA 1 ^ 3 fm, in accord with the liquid drop 
model. Each nucleon occupies a volume of h 3 , so that phasespace is uniformly filled. The 
initial momenta are randomly choosen between and the local Thomas-Fermi-momentum. 
The Ap and At nucleons interact via two- and three- body skyrme forces, a Yukawa potential 
and momentum dependent interactions. Subsequently, the FMD |ft8 |, AMD |39[ and PQMD 



models have been developed to offer an improved treatment of the Pauli-principle. 
Isospin is treated explicitely (in the socalled "I" QMD version), a symmetry potential (to 
achieve corrected distributions of protons and neutrons in the nucleus) and explicit Coulomb 
forces between the Zp and Zt protons are included. 

Pion production is treated via the delta resonance P^TTpS]. A frozen A approximation 



(infinite lifetime for the A resonance) had been used in other versions of the QMD model. 

The hadrons are propagating under the influence of the potential in Hamilton's equations 
of motion: 
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Pi = 



dH 

dqi 




dH 



with 



H = T + V = 



2rrii 



+ fiiff V, t) V l] fj{f\ p', t) dfdf'dp dp 



i j>i 



The baryon-potential V™ = V& + V?uk + VcU + + V tlm consists of 



v ioc = hS(ri - fj) + fk(r'',p'')t 2 5(fi-f j )5(f i -f k )df''dp' 



k>j>i 




V sym = t 6 —T 3i T 3j 6(fi - fj) 



The three-body term in V^ c as stated above is only valid for a hard equation of state (for a 
soft equation of state a VUU-type formulation ~ p 7 has to be used), Zi, Zj denote the charges 
of the baryons % and j and T 3i , T 3 j are their respective T 3 components. The meson-potential 
only consists of the Coulomb potential. 

The parameters /i and t 1 ...t 6 are adjusted to the real part of the nucleon optical po- 
tential. For the density dependence of the nucleon optical potential standard Skyrme type 
parametrizations are used. Two different equations of state have been implemented: A hard 
equation of state (H) with a compressibility of 380 MeV and a soft equation of state (S) with 
a compressibility of 200 MeV P~S|J5T[ ] . A fit of the momentum dependence to measurements 
f4~3| , |4"4|1 of the real part of the nucleon optical potential H6j |2~l],[l5| yields: 



v ' \PoJ 

The equation of state in its standard Skyrme type parametrization including momentum 
dependence then reads: 




The mean field notation with the parameters a, /3, 7, 8 and e has been chosen for rea- 
sons of simplicity and in order to compare the parameters with those used in VUU/BUU 
calculations. Their values can be found in table |I[ While the forces are calculated via the 
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nucleoli density in VUU/BUU calculations, a sum over two-particle interactions is performed 
in QMD/IQMD calculations. 

The parameters ^...tg are calculated in the IQMD model before the initialization of the 
projectile and target nuclei from the tabulated values of a, (3, 7, 5 and e which serve as input. 
The width and normalization of the Gaussian wave-packets have to be taken into account 
for the proper determination of the force-parameters. 

Hard N-N-collisions are included by employing the collision term of the well known 
VUU/BUU equation . The collisions are done stochastically, in a similar way 



as in the cascade models [48 , 49| . Two particles collide if their minimum distance d fulfills 



d < d = J—, otot = v(y/s, type). 

V 7T 

"type" denotes the ingoing collision partners (N — N, N — A, N—tt...). In addition, the Pauli 
blocking (of the final state) of baryons is taken into account by checking the phase space 
densities in the final states of a two body collision. The final phase space fractions Pi and P2 
which are already occupied by other nucleons are determined for each of the two scattering 
baryons. The particular attempt for a collision is then blocked with the probability 

Puock = 1-(1-Pa)(l-P 2 ) 

Whenever an attempted collision is blocked the scattering partners maintain the original 
momenta prior to scattering. Delta decays are checked in an analogous fashion with respect 
to the phase space of the resulting nucleon. 

Pions are formed in the IQMD model via the decay of the delta resonance. The following 
inelastic reactions are explicitly taken into account and constitute the imaginary part of the 



pion optical potential, which is dominant in the 1 GeV/u energy domain [50 

a) N N — > AN (.hard-delta-production) 

b) A -> Nit (A-decay) 

c) AN -> NN (A-absorption) 

d) Nn — > A (soft-deita-production) 

Elastic 7T — 7r,7r — N, n — A, A — A, A — N scattering is not taken into account. Experimental 
cross sections are used for processes a) and d) pi] , as well as for the elastic N-N-collisions. 
The respective cross sections are shown in figure [I]. 



For the delta absorption, process c), we use a modified detailed balance formula [52 



The conventional detailed balance formula is only correct for particles with infinite life- 
times (zero width). If the principle of detailed balance is applied to the delta resonance, 
then its finite width has to be taken into account: 
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da NA^NN ± da NN^NA 



dfl 4 pata dQ 

with 



— I p NA A r (M)2MdM 



A r (M) - 1/2 



(M-m A ) 2 + (r/2) 2 
The mass-dependent A-decay width has been taken from ||53||: 



Fa(p-k) = ^r r o with r(p) - 



r(po) 1+M + (-£■ 



is the decay momentum of the pion, p = 227 MeV/c, p\ = 238 MeV/c, p 2 = 318 MeV/c 
and T = 120 MeV. The A decays isotropically in its restframe. 

The elastic nucleon-nucleon scattering angular distribution is taken to be [54| : 

where t is — q 2 , the squared momentum transfer and 

(3.65 (y/s- 1.8766)) 6 



A(s) = 6 



1 + (3.65 (y/s- 1.8766)) 6 

y/s is the cm. energy in GeV and A is given in (GeV/c)~ 2 . 

The inealstic channel is treated in an analogous fashion. The parametrization suggested 



by Huber and Aichelin |55| is used: fitted differential cross sections are extracted from OBE 
calculations: 

j* 1 ~ a(s) exp(b(s) cos 60 , 
ail 

a(s) and b(s) are functions of y/s and vary in their definition for different intervals of y/s 
(see table [H]). 

Pions propagate between collisions (imaginary part of the pion optical potential) on 
curved trajectories with Coulomb forces acting upon them. The different isospin channels 
are taken into account using the respective Clebsch-Gordan-coefficients: 

A++ -> 1(p + tt+) A+ -> f(p + 7r°) + i(n + 7r+) 

A -> |(n + 7r°) + i(p + 7r-) A- -> l(n + 7r-) 

The real part of the pion optical potential is treated in the following manner: As far as 
the pion is bound in a A-resonance, the density and momentum dependent real part of the 
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nucleoli optical potential is applied as an approximation to the (yet unknown) real part 
of the A optical potential. Due to the large tt — N cross section, intermediate pions are 
quite frequently bound in a Delta resonance and in that intervals the real part of the pion 
optical potential is substituted by the real part of the A optical potential. Free intermediate 
and final charged pions experience coulomb forces which contribute to the real part of 
the pion optical potential. Recent investigations on the influence of the nuclear medium 
correction to the dispersion relation of the free pion have shown conflicting results p>6]J57 



with respect to the importance of the modification for low momentum pions. However, both 
calculations show that the high energy part of the pion spectrum remains unchanged by this 
modification. Since our results are mainly for this high energy contribution we omit this 
medium correction, until a consensus has been achieved on the proper form of the respective 
medium contribution. 

After a pion is produced (be it free or bound in a delta), it's fate is governed by two 
distinct processes: 

1. absorption tc N N ^ AN ^ N N 

2. scattering (resorption) it N — > A — ► n N 

In the CASCADE mode all real forces are turned off: nucleons, pions and deltas are propa- 
gated on straight lines between collisions. 

III. FORMATION AND SPECTROSCOPY OF A RESONANCE MATTER 



Recently an old subtopic |T0| , [58[| of this research has received renewed attention |59H63|: 



The possibility of producing A-matter (or in more general terms: resonance matter). At 
beam energies above a few hundred MeV/nucleon, the nucleons can be excited into A- 
resonances. If the density of these resonances is as high as the nuclear matter ground state 
density, then a new state of matter, A-matter, has been created. One of the potential signals 
for the presence of A-matter is the creation of pions as decay-products of the A-resonance. 

How can A-matter be produced? Figure |2| shows the pion - nucleon cycle in the IQMD 
model. The scheme describes (for impact parameters b < 5 fm and averaged over 60 fm/c 
possible processes linked to the creation of A-matter. The probabilities in the boxes always 
refer to the vertices they are directly connected with. A-resonances are initially produced 
via inelastic nucleon nucleon scattering. The produced resonances can either be reabsorbed 
via inelastic scattering or decay by emitting a pion. The pion can then either freeze out or 
interact with a nucleon to form a A again. In case the A has been absorbed the corresponding 
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high energetic ergetic nucleon might have a second chance of becoming a A by inelastic 
scattering. It could also transfer energy via elastic scattering onto another nucleon which 
then could scatter inelastically and form a new A. A nucleon interacts in the average about 
three times before it freezes out. This value fluctuates considerably, depending on whether 
the nucleon was in the participant zone (geometrical overlap of the colliding heavy ions) or 
in the spectator zone of the collision. 

Unfortunately, the probablity for a nucleon to undergo inelastic scattering and to form 
a A during the heavy ion collision is as low as 10%. The main process for sustaining A- 
matter is the A — > Nit —* A loop, which, however, first has to be fueled by the NN — > AN 
process. The average pion passes approximately three times through this loop (it has been 
created by the decay of a hard A). However, 30% pass more than 6 times through the loop. 
For nucleons the probability of forming a soft A i.e. via ttN — > A is almost twice as high 
(A-matter pump) than the probability of forming a hard A via NN — > NA. 

Figure ^|a) shows the time evolution of the total baryon, nucleon and A densities in units 
of p/po (top). The densities are calculated in a sphere of 2 fm radius around the collision 
center. Between 5 fm/c and 20 fm/c more than 20 A-resonances can be found in the system: 
This time interval coincides with the hot and dense reaction phase. At 10 fm/c up to 55 
resonances are present in the total reaction volume (keep in mind this is not in the 2 fm 
test sphere). A A multiplicity of > 40 can be sustained for an interval of 10 fm/c, 6 times 
longer than the lifetime of a free A-resonance. However, this is not pure A-matter. in the 
small test volume shown in figure 2a the resonance density is 0.5 po and the nucleon density 
is 2.2 pq: the A-contribution is 20% in the test volume which contains, as a matter of fact, 
only 2.5 resonances. The total multiplicity of A resonances is just about 10% of the total 
nucleon multiplicity. 

However, it is obvious that the other As can be distributed all over the reaction volume. 
Figure |||b) shows the A density distribution as experienced by the A's in the system at 5, 
10 and 20 fm/c. The densities were calculated by summing over all contributing Gaussians 
of all As in the system at the locations of the respective As. We would like to point out 
that the mean A-density experienced by the As is about 0.25 po- Less than 1% of the As 
experience A densities around 0.5 po- However, enough As are in the system to show signs 
of collectivity such as collective flow in the reaction plane. Its measurable signiture (the 
pion p* x {y) distribution in central collisions) will be discussed in one of the follwing sections. 



S 



IV. INCLUSIVE PION OBSERVABLES 



This section deals with inclusive pion production in relativistic heavy ion collisions. 
Figure |] shows the predicted impact parameter dependence of the multiplicity of n~ , 7r° and 
7r + for Au+Au collisions at 1 GeV/nucl. incident beam energy. A hard EoS with momentum 
dependent interaction is used in the calculation. For central collisions (6 = fm) the total 
pion multiplicity is approximately 55. However, the average pion multiplicity is about 19 (8 
7r~, 6 7T° and 5 rr + ) for a minimum bias impact parameter distribution. 

The mass dependence of the total pion multiplicity is shown in figure For light collision 
systems the multiplicity increases linearly with the system mass. However, for heavier 
systems the increase is less than the linear extrapolation, this is due to pion absorption. 
The values are for b = fm calculations of the systems 20 Ne+ 20 Ne, 40 Ca+ 40 Ca, 58 Ni+ 58 Ni, 

93 Nb+ 93 Nb and 197 Au+ 197 Au> 

The polar angular distribution — for tt~ , tt° and tt + is shown in Figure |6| for mini- 
mum bias events (a) and for tt~ in central vs. minimum bias events in (b). A horizontally 
flat distribution would correspond to isotropic emission. For minimum bias events (top) a 
strong peaking towards forward-backward angles is observed, most prominently for tt~. 

This is important for the extrapolation of total yields from spectra measured at $ c .m. = 
90° - if the midrapidity spectra are used to extrapolate (with the assumption of a flat dis- 
tribution) the total yield may be underestimated by a factor of 2. The anisotropy decreases 
when studying central collisions (bottom). This dependence can be explained by the decay 
of A-resonances in the projectile- and target- spectator regions. The difference betweeen 
the distributions of ir~ and ir + also results in a forward-backward peaking of the iV 7r -/iV 7r + 
ratio. This phenomenon has already been experimentally observed for light collision systems 



at the BEVALAC 64 



Figure [7] shows a comparison of inclusive tt° spectra for Au+Au and Ca+Ca (minimum 
bias and y c . m . ± 0.16) between the IQMD model and data published by the TAPS collabo- 
ration ||65|| . Whereas the model shows reasonable agreement with the heavy system Au+Au 
it overpredicts the 7r° yield of the light system Ca+Ca by approximately 60%. The charged 
pions (for Au+Au collisions) are shown in figure |8| together with ir + data from the KaoS 
collaboration [B6|. The slope of the 7r + spectrum in the model calculation agrees well with 
the KaoS measurements. However, the multiplicity as predicted by the model is approxi- 
mately 20% above the KaoS measurements. Both, the calculation and the measurements, 
have been acceptance corrected to the rapidity interval y c . m . ± 0.16 and may directly be 
compared to Figure [7|. Especially the measurement of high energy pions is of great interest. 
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They correlate directly to early freeze out times and heavy A resonances |67] . 

The mass dependence of pion production and its sensitivity towards the transverse mo- 
mentum p t can be studied more clearly by plotting the ratio of the p t spectra for Au+Au 
and Ca+Ca versus the transverse momentum p t (Figure ||]). A comparison between data 
from the TAPS collaboration |6f| and the IQMD model is shown in figure [|. The model 



underpredicts this ratio by approximately a factor of 2, In particular for low transverse mo- 
menta . However, this holds for all transport model calculations which have been compared 



to the TAPS data [65]. A previous comparison between the IQMD model and the TAPS 
data in ref. []65| showed far larger disagreement between model and data. This was due to 
an improper normalization of the theoretical calculations. 

The yield of low p t pions in the heavy system is underpredicted by 10%. This is a 
common problem of most transport theories dealing with heavy ion collisions. Suggested 
explanations for this underprediction include in-medium effects of pions in nuclear matter 
|57j and the neglect of Bose-enhancement due to the bosonic nature of the pions. 

The importance of the inclusion of Coulomb forces and energy dependent tt — N cross 
sections can be shown by plotting the tt~ to ir + ratio versus the transverse momentum p t 



(Figure [H|). The solid line shows the full calculation including Coulomb forces. For high p t 
the n~ /n + ratio decreases towards 1, whereas for low pt it increases to 2.5 - considerably 
higher than the value of 1.8 predicted by the A-isobar model. The dashed line shows 
a calculation without Coulomb forces. This ratio remains constant at 1.8. The (small) 
remaining variations might be due to the different energy dependence of the tt + — p and 
7r~ — p inelastic cross sections. 

V. PION NUCLEON CORRELATION IN THE EVENT-PLANE 

The hydrodynamical model predicts a bounce-off of nuclear matter in the reaction plane 
which has experimentally indeed been discovered [5"9,70"|. The bounce-off is depicted 



by plotting the in-plane transverse momentum p~* x (y) versus the rapidity y. For nucleons 
and light fragments a horizontal s-shape is typically seen with negative p x {y) values for 
y < y c .m. and positive p" x (y) values for y > y c , m ,. Figure [11] shows the p~x(y) distribution 



for 7r + and protons in Au(lAGeV)Au collisions with a minimum bias impact parameter 
distribution. The protons show the expected collective flow []7l]-[[3|]. The p x of the pions, 
however, is anticorrelated to that of the protons. A similar proton - tt anticcorelation has 
been measured for the asymmetric system Ne+Pb at 800 MeV/nucl. by the DIOGENE 
collaboration Transport model comparisons to the DIOGENE data with the IQMD 
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24| and the BUU |25| model have shown good agreement with the data. 

We have studied the origin of the particular shape of the pion angular distribution and 
the p x spectrum by sequentially suppressing first the sq/£-<ieZta-production and then the 
delta-absorption (while allowing the so/t-de/ta-production). If we deactivate the soft-delta- 
production (see Figure 3), tt N — > A, pions are neither scattered nor absorbed after the 
initial production. No p x for pions is observed. In order to decide whether the p x spectrum 
is caused by absorption or by scattering we now deactivate the reaction AN — > NN. 
We thus suppress pion absorption but allow scattering - the anticorrelation between pions 
and protons in the p" x returns. In contrast to previous publications, which investigated 
the asymmetric system Ne(800AMeV)Pb and suggested the anticorrelation of pionic and 
nucleonic p" x at target rapidities to be caused by pion absorption |]2"5| , our investigation 
reveals the p" x spectrum of the pions to be dominated by the pion scattering process [|6|]. 

The following simplified picture can explain the origin of the observed phase space dis- 
tribution: The A decays isotropically in its rest-frame, therefore 50 % of the pions are 
emitted with a positive p x and 50 % with a negative p x . At target rapidity those pions 
which obtain a positive p^-value usually do not have the chance to rescatter: Most of the 
target nucleons are located in the negative p x area! Those pions which do rescatter at target 
rapidity are the ones with an initially negative p x : Every time a A decays (isotropically) 
there is a 50% chance that this pion is emitted upward, i.e. into an azimuthal angle between 
—90° < < 90°. These 0-values characterize the hemisphere of positive p x , by definition. 
This leads for^ 50% of the pions with - originally - negative p x to a shift towards a positive 
p x . This remains true even after transforming back into the laboratory frame. The same 
consideration applies vice versa for projectile rapidity: Most projectile nucleons are located 
in the positive p x area. The pions are rescattered in this area which results in a negative p" x 
and a maximum in the azimuthal angular distribution in the 90° < <p < 270° interval. 

Figure shows the in-plane transverse momentum p" x versus rapidity y (in the cm. 
system) for 7r + in central collisions of Au+Au (with impact parameters b< 3 fm). In contrast 
to (semi-) peripheral collisions, however, p" x is correlated for pions and nucleons in central 



collisions because of the bounce-off of A- resonances | 75| . The square markers in Figure 
|i"2] depict a calculation with the hard equation of state without momentum dependence, 
the circles show the same equation of state including momentum dependence whereas the 
triangles represent a CASCADE calculation, i.e. a non-equilibrium free gas. 

The momentum dependence enhances the p" x of the pions. This effect is due to the 
bounce-off of the A resonances J7^] which in our model is enhanced because the momentum 
dependence for the A resonances is included in the same way as for the nucleons. 
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The CASCADE calculation, however, shows the opposite behaviour. The p" x of the pions 
has negative sign to that of the calculations with the density dependent equations of state. 
This behaviour can be explained by the lack of hadron collective flow in CASCADE cal- 



culations f76|. The pions would then be expected to be emitted isotropically (px(Y) = 0). 
However, pion scattering from small caps of spectator matter being present at impact pa- 



rameters around 3 fm causes the observed anticorrelation ||26|| . In order to investigate the 
density dependence of the nuclear equation of state and in order to show the differences be- 
tween CASCADE calculations and calculations including the equation of state more clearly 
we use the robust observable pf r which for nucleons is defined as 

A P +A T 

E Pl ■ sgn(?/i - y c . m ) 

dir 8=1 



Px 



A T + A P 

(the adaptation for pions is straightforward) and plot it versus the impact parameter (Figure 
|13D . For positive values of p d " the pion p" x vs. rapidity distribution is correlated to that of 
the nucleons. For negative values an anticorrelation is observed. 

Figure shows the respective calculations for the hard and soft equations of state 
(including momentum dependence) and for the CASCADE calculation. For small impact 
parameters the calculations with equation of state show a correlation between pion and 
nucleon bounce-off. At semiperipheral impact parameters we observe a sign reversal. As 
mentioned above, the anticorrelation between nucleon and pion bounce-off is caused by 
pion scattering in spectator matter ||26|| . In contrast, the CASCADE calculation exhibits 
a negative p d x tr for the whole impact parameter range. The momentum transfer pf r shows 
a systematic difference between the hard and soft equation of state. However, very high 
statistics and high precision impact parameter classification are necessary to experimentally 
exploit this sensitivity towards the determination of the nuclear equation of state. The 



results of figures |Tl] and [L3| show clearly that even in the domain of particle production 
(ir, K, T), p, p, uj) CASCADE simulations predict distinctly different phase space distributions 
for baryons and mesons at central impact parameters. 



VI. AZIMUTHAL CORRELATIONS OF PIONS PEPENDICULAR TO THE 

EVENT PLANE 

Now let us investigate particle emission perpendicular to the reaction plane. The hy- 
drodynamical model predicted a squeeze-out of high energetic nucleons perpendicular to 
the reaction plane [P,|i|,|77|. This effect, which has also been predicted by QMD-calculations 
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2^ , [78| -pT| and has been confirmed by experiment ||B2] -|84| , is due to the high compression of 



nuclear matter in the central hot and dense reaction zone (it is a genuinly collective effect, 
increasing linearly with A). 

Do pions show a similar behaviour? The azimuthal (tp) distribution of the pions is plotted 
to investigate this question, (p is the angle between the transverse momentum vector p t and 
the x-axis (which lies in the reaction plane and is perpendicular to the beam axis). Thus 
if = 0° denotes the projectile hemisphere and if = 180° corresponds to the target hemisphere. 



Figure [bj shows the respective distributions for neutral pions in the transverse momen- 
tum bins pt < 50 MeV and pt > 400 MeV at a minimum bias impact parameter distribution. 
The distributions have been normalized in order to fit into the same figure. The analysis was 
performed from 0° to 180° and then symmetrized for 180° to 360°. The plotted distributions 
have been extracted by fitting the calculated points (shown for the high p t bin) according to 
the function a ■ (1 + b ■ cos(0) + c ■ cos(20)). The azimuthal angular distribution for 7r° with 
low pt shows maxima at (p = 0° and tp = 180° corresponding to a preferential emission in the 
reaction plane. The high p t it , however, show a maximum at <p = 90°. This maximum is 
associated with preferential particle emission perpendicular to the reaction plane. The inlay 
shows data from the TAPS collaboration for the region 400 MeV < pt < 600 MeV and 
midrapidity. We observe a good qualitative agreement between the theoretical prediction 
and the experiment. It should be noted, however, that both, theory and experiment, need 
much better statistics to allow a conclusive quantitative comparison. 

The magnitude of the observed anisotropy and its dependence on impact parameter and 
transverse momentum is best studied by using the following ratio: 



R, 



g (^ = 90°) + f (^ = 270° 



out/in dN , _ „ - dN 



y=yc 



For positive R ou t/in values pions are emitted preferentially perpendicular to the reaction 
plane. Figure |T3| shows the transverse momentum dependence of R ou t/in for Au+Au collisions 
with an impact parameters bewteen b=5 fm and b=10 fm and at midrapidity (y c . m . = ±0.2): 
In contrast to pions with low transverse momentum, which are emitted preferentially in 
the reaction plane, high p t pions are preferentially emitted perpendicular to the reaction 
plane. This effect is stronger for 7r + than for 7r~. The difference is due to the different 
ttN — > A production cross section for 7r + and n~ and due to Coulomb forces pushing the 7r + 
away from the spectator matter which is located mostly in the reaction plane. The 7r~ on 
the other hand are being attracted by those spectator-protons. These effects decrease the 
number of tt~ leaving the reaction zone in a direction perpendicular to the reaction-plane. 



13 



However, the statistics accumulated so far are not large enough for a more detailed study 



of the differences between positive and negative pions. The inlay of Figure [15] shows recent 
measurements from the KaoS collaboration which confirm the predicted systematics of 
the pt dependence. Imposing the limited acceptance of the KaoS spectrometer on the IQMD 
calculations would reduce the available statistics by one order of magnitude. Experimental 
uncertainties in the determination of the proper reaction plane result in a reduction of the 
measured Rout /in values which are difficult to compensate. Therefore a direct quantitative 
comparison between the KaoS measurements and IQMD calculations is not feasable at this 
point of time. 

We have investigated the cause of the observed preferential emission perpendicular to 
the reaction plane: Pion absorption as well as scattering can be eliminated by deactivating 
the reaction n N — > A, then no sqeeze-out is observed. 

In order to decide whether the anisotropy is caused by absorption or by scattering the 
reaction A N N N can be deactivated. Thus pion absorption is suppressed but scat- 
tering is allowed: no anisotropy is observed. Therefore we conclude that the anisotropy is 
dominated by the pion absorption process 



Figure [16] shows the distribution of the number of delta generations wa a pion goes 
through before its freeze out. Here wa is shown for tt + emitted both in the reaction plane 
as well as perpendicular to it. (wa — 1) is therefore the number of times a pion scatters 
before freeze out. We observe that 90% of the produced pions scatter at least once before 
leaving the reaction zone. A large number of pions scatters even more often, 2% up to 10 
times! The observed preferential emission perpendicular to the reaction plane is due to an 
excess of high p t pions which on the average have undergone fewer collisions (< 2) than the 
pions in plane. Those pions which make this effect do rescatter rarely, they are emitted early 
but carry information on the high density phase of the reaction. They stem from the decay 
of the most massive delta resonances which are mostly produced early on in the reaction 
|67|. Therefore high p t pions emitted perpendicular to the event-plane should be the most 



sensitive pionic probes for the investigation of the hot and early reaction zone. 

Figure [T7| shows the azimuthal angular distribution of high p t (pt > 400 MeV) neutral 
pions at midrapidity and impact parameter b=6 fm. The different curves show calculations 
for hard (circles) and soft (squares) equations of state (including momentum dependence) 
and a CASCADE calculation (triangles), (p is the angle between the transverse momentum 
vector p t and the x-axis (which lies in the reaction plane and is perpendicular to the beam 
axis). The out-of-plane pion squeeze-out is clearly seen by the pronounced maximum at 
<f = 90° for both equations of state. To enhance the statistics all particles are projected 
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into the 0° < ip < 180° hemisphere. The full and dashed lines are least square fits with the 
function f{ip) = a(l + siCOs(<p) + s 2 cos(2y9)) which has been used to fit the squeeze-out 
phenomenon |8"!2| , |8"3|] . The curves show an extrapolation to the full azimuthal angular range. 
The distributions are normalized per particle in order subtract the influence of different 
equations of state on the pion multiplicity. Within error-bars both equations of state exhibit 
the same out-of-plane pion squeeze-out. There is a trend for the hard equation of state to 
exhibit an enhanced out-of-plane pion squeeze-out but this trend might be too small to be 
useful for an experimental distinction between the different equations of state, whereas - in 
contrast - the in-plane pion bounce-off shows a clear difference for the two cases (see above). 
The CASCADE calculation does not exhibit any significant out-of-plane pion squeeze-out 
for b=6 fm! However, for larger impact parameters also CASCADE calculations exhibits a 
pronounced out-of-plane pion squeeze-out (Figure fL8|) . 

Figure [18] shows R ou t/in versus impact parameter b for the hard equation of state with 
and without momentum dependence of the real part of the nucleon optical potential. The 
momentum dependence causes a drastic increase of R ou t/in for impact parameters larger than 
3 fm. In the CASCADE calculation the onset of the out-of-plane pion squeeze-out is shifted 
toward larger impact parameters as compared to the calculations including the equation of 
state (see also Figure |T7|) . For peripheral collisions R ou t/in reaches the same magnitude for 
the CASCADE calculation and the hard equation of state without momentum dependence. 



VII. CONCLUSION AND OUTLOOK 



The physics responsible for the in-plane pion bounce-off (pion scattering) and the out- 
of-plane pion squeeze-out (pion absorpion) differs completely from the compressional effects 
governing the in-plane nucleon bounce-off and out-of-plane nucleon squeeze-out. It is the 
pion-nucleon interaction which creates the sensitivity towards the density and momentum 
dependence of the nucleon optical potential. Therefore it is understandable that we observe 
a strong sensitivity towards the equation of state in the reaction plane whereas the sensitivity 
towards the equation of state perpendicular to the reaction plane is limited to the momentum 
dependence: The (anti-) correlation in-plane is caused by multiple pion nucleon scattering 
PS| with the bounced-off nucleons, which show a strong sensitvity towards momentum 
and density dependence IftlH . The pion squeeze-out perpendicular to the reaction plane, 
however, is dominated by high p t pions which have undergone less rescattering than those in 
the reaction plane ||28|| . The abundance of these high p t pions is correlated to the multiplicity 
of high p t nucleons which increases if the density dependence is included. 
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We have investigated the dependence of pionic in-plane bounce-off and out-of-plane 
squeeze-out on the nuclear equation of state. A strong sensitivity towards the density 
dependence is observed for the in-plane pion bounce-off whereas the out-of-plane pion 
squeeze-out shows only a small sensitivity. Both effects show a strong sensitivity toward the 
momentum dependence. CASCADE calculations - which we see as a crude approximation 
to QMD - give different phase space distributions for pions in both cases. It should be easy 
to resolve experimentally these two clearly qualitatively different distinct scenarios. The 
determination of the equation of state will require - on the other hand - a more sensitive 
(and sensible) quantitative comparison to theory, including a improved treatment of the A 
and pion optical potential. 

The nuclear equation of state cannot be extracted from one observable alone. All ob- 
servables known to be sensitive to the equation of state have to be fitted simultaneously by 
the respective model in order to claim success. In this article we have added additional, here 
pionic, observables which have to be taken into account for obtaining the final goal: The 
nuclear equation of state. 
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TABLES 



a (MeV) 


(3 (MeV) 


7 


6 < MeV ) £ (oSf 


S -356 
SM -390 
H -124 
HM -130 


303 
320 
71 
59 


1.17 
1.14 
2.00 
2.09 


1.57 500 
1.57 500 


TABLE I. Parametersets for the nuclear equation of state used in the IQMD model. S and 
H refer to the soft and hard equations of state, M refers to the inclusion of momentum dependent 
interaction. 


x = v« (GeV) 


a (fm) 




b 


2.104 - 2.12 
2.12 - 2.43 
2.43 - 4.50 


294.6 (x - 2.014) 2 - 578 

0.01224 
(x-2.225) 2 +0. 004112 

(2.343/x) 43 - 17 




19.71 (x - 2.014) 1 - 551 
19.71 (x - 2.014) 1 - 551 
33.41 arctan(0.5404 (x - 2.146) - 9784 ) 



TABLE II. a(s) and b(s) as functions of the cm. energy. 
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FIGURES 
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3.5 

s ' / GeV 

FIG. 1. Tabulated cross sections in the IQMD model. The upper frame shows the total, elastic 
and inelastic proton neutron cross section. The middle frame compares the total proton neutron 
with the total proton proton cross section and the lower frame shows the total 7r + proton (or 7r~ 
neutron) cross section. The other pion nucleon cross sections are determined by scaling the ir + 
proton cross section either with 1/3 (tt + neutron and 7r~ proton) or with 2/3 (ir° neutron and 7T° 
proton). 
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FIG. 2. Pion - nucleon cycle in the IQMD model. The scheme describes (for b < 5 fm 
and time-averaged) all in the model possible processes linked to the creation of A-matter. The 
probabilities in the boxes always refer t o the vertices they are directly connected with. The main 
process for sustaining A-matter is the A — ► Ntt — ► A loop, which, however, first has to be fueled 
by the NN -> AN process. 
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FIG. 3. Time evolution of the the total baryon, nucleon and A-resonance density in units 
of p/po (a) and A density distribution the respective As experience for 5, 10 and 20 fm/c. The 
densities in the upper frame (a) are calculated in a sphere of 2 fm radius around the collision center. 
The hot and dense reaction phase lies between 5 and 20 fm/c during which approximately 10% of 
the nucleons are excited to A-resonances. Up to 50% nucl. matter ground state density is reached 
by the A-resonances. However, less than 1% of the As experience such high densities. The average 
A density which is felt by the As is approximately 0.25 po at 10 fm/c. The densities in the lower 
frame (b) were calculated by summing over all contributing Gaussians of all A's in the system at 
the locations of the respective A's. 
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impact parameter b (fm) 
FIG. 4. Multiplicity of ir~,ir° and ir + versus impact parameter b for Au+Au collisions at 1 
GeV/nucl. incident beam energy. A hard EoS with momentum dependent interaction is used. For 
central collisions (6 = fm) the total pion multiplicity is approximately 55. For a minimum bias 
impact parameter distribution the average pion multiplicity is about 19 (8 ir~ , 6 ir° and 5 ir + ). 
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FIG. 5. Total pion multiplicity versus collision system mass at 1 GeV/nucl. beam energy. For 
light collision systems the multiplicity increases linearly with the system mass. However, for heavy 
systems the increase is less than linear due to pion absorption. The values plotted were extracted 
from b = fm calculations of the systems Ne+Ne, Ca+Ca, Ni+Ni, Nb+Nb and Au+Au. 



27 



0.12 




« 0.035 

• r—i 

2 



o3 



0.03 



0.025 
aj 0.02 



fH 

o 
d 



0.015 



§ 0.01 

CO 

g 0.005 



0.0 




7r minimum bias 
7r central, b< 3 fm 



b 



45 



0, 



90 



135 



180 



FIG. 6. Polar angular distribution -j — ^¥ — for 7r , 7r° and 7r + in minimum bias (a) and for 

° d COSt) cm. v ' 

pi~ in minimum bias and central (b) Au+Au collisions. A horizontally flat distribution would 
correspond to an isotropic emission. For minimum bias events (top) a strong peaking towards 
forward-backward angles is observed, most prominently for 7r~. The anisotropy decreases when 
studying central collisions (bottom). It can be explained by the decay of A-resonances in the 
projectile- and target- spectator regions. 
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FIG. 7. Comparison of inclusive 7r° spectra for Au+Au and Ca+Ca (minimum bias) 
collisions between the IQMD model and data measured by the TAPS collaboration. A hard EoS 
including momentum dependence is used and a rapidity cut according to the acceptance of the 
TAPS spectrometer is employed. The model shows reasonable agreement with the heavy system 
Au+Au but it overpredicts the 7r° yield of the light system Ca+Ca by approximately 60%. The 
yield of low p t pions in the heavy system is underpredicted by 10%. 

FIG. 8. Inclusive n" and ir + spectra for Au+Au (minimum bias) collisions at 1 GeV/nucl. 
as calculated with the IQMD model and a ir + spectrum measured by the KaoS collaboration. A 
hard EoS including momentum dependence is used in the IQMD calculation and a rapidity cut 
according to the acceptance of the TAPS spectrometer is employed; the KaoS measurement has 
been acceptance-corrected. The slope of the ir + spectrum in the model calculation agrees well with 
the KaoS measurements. However, the multiplicity as predicted by the model is approximately 
20% above the KaoS measurements. 

FIG. 9. Ratio of the pion yield from Au+Au and Ca+Ca collisions plotted versus the trans- 
verse momentum p t . The figure shows a comparison between data measured by the TAPS col- 
laboration and an IQMD calculation (minimum bias, hard equation of state with momentum 
dependence). For low transvere momenta the model underpredicts the data approximately by a 
factor of 2. 

FIG. 10. 7r~ to 7r + ratio versus transverse momentum p t for Au+Au collisions (minimum 
bias) at 1 GeV/nucl. with a hard EoS and momentum dependence. The solid line shows the 
full calculation including Coulomb forces. For high pt the ratio decreases towards 1 whereas for 
low pt it increases to 2.5 - considerably higher than the value of 1.8 which the A-isobar model 
would suggest. The dashed line shows a calculation without Coulomb forces acting upon the pions. 
Within the errorbars the ratio remains constant at a value around 1.8. 

FIG. 11. Rapidity y vs. p" x /m for tt + and protons in Au+Au collisions at 1 GeV/nucl. with 
minimum bias impact parameter distribution and a hard EoS including momentum dependence. 
The protons show the expected bounce-off. The p" x of the pions is directed oppositely to that of 
the protons. This effect is caused by rescattering of pions from large chunks of spectator matter. 



29 



FIG. 12. in-plane transverse momentum p* x versus rapidity y (in the cm. system) for ir + in 
central Au+Au collisions with impact parameters b< 3 fm. Calculations with a hard equation of 
state without momentum dependence (squares), the same equation of state with momentum de- 
pendence (circles) and a CASCADE calculation (triangles) are shown. The effect of the momentum 
dependence is considerable, exhibiting the sensitivity of p x on the baryon flow. The CASCADE 
calculation gives a different phase space distribution due to its lack of collective baryon flow. 

FIG. 13. p^ r versus impact parameter b for Au+Au collisions at 1 GeV/nucl. with hard 
and soft equations of state, both with momentum dependence, and for the CASCADE calculation. 
Note the clear sensitivity for the equation of state on p^ r . The CASCADE calculation exhibits 
an anticorrelation between pions and nucleons for the whole impact parameter range, due to pion 
nucleon scattering and its lack of collective baryon flow. 

FIG. 14. Normalized azimuthal angular distribution dN/dip for ir° with low and high trans- 
verse momentum p t at mid-rapidity in the reaction Au(lAGeV)Au with minimum bias impact pa- 
rameter distribution, a hard equation of state and momentum dependent interaction. The points 
were fitted according to the function a ■ (1 + b ■ cos(</>) + c • cos(20)). The maximum at </? = 90° 
corresponds to a preferential emission of high pt pions perpendicular to the reaction plane. This is 
due to pion absorption by large pieces of baryonic spectator matter located predominantely in the 
reaction plane. Perpendicular to the plane there is no such spectator matter and pions with high pt 
can leave the reaction zone without further interaction. Low pt pions have rescattered more often 
which is only possible in the reaction plane. The inlay shows data from the TAPS collaboration 
for the region 400 MeV < pt < 600 MeV and midrapidity. 

FIG. 15. Squeeze-out ratio R ou t/in versus transverse momentum for ir + and 7r _ . Pions with 
Pt > 200 MeV are preferentially emitted perpendicular to the reaction plane. Pions with pt < 100 
MeV are emitted isotropically because they have undergone frequent rescattering which can only 
happen due to spectator matter in the reaction plane. The differences between tt~ and n + are due 
to Coulomb forces. The inlay shows data on ir + from the KaoS collaboration. 

FIG. 16. Distribution of the number of delta generations i%a a pion goes through before its 
freeze out for ir + emitted in the reaction plane and perpendicular to it. 90% of the produced 
pions scatter at least once before leaving the reaction zone. The observed preferential emission 
perpendicular to the reaction plane is due to an excess of pions which on the average have undergone 
fewer collisions (< 2) than the pions in plane. 
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FIG. 17. Azimuthal angular distribution dN/dip for neutral pions calculated with hard and 
soft equations of state (both with momentum dependence) and a CASCADE calculation. All 
calculations were performed with an impact parameter of b=6 fm. Both equations of state exhibit 
approximately the same angular distribution whereas the CASCADE calculation does not exhibit 
any peak perpendicular to the event plane. For larger impact parameters, however, also the 
CASCADE calculation shows a pronounced squeeze-out. 

FIG. 18. squeeze-out ratio R ou t/in versus impact parameter for neutral pions. The calculations 
were performed with the hard equation of state with (circles) and without (squares) momentum de- 
pendence as well as in CASCADE mode (triangles). Cuts around mid-rapidity (—0.2 < y c . m . < 0.2) 
and for high transverse momentum (p t > 300 MeV) were employed. For large impact parameters 
the CASCADE calculation agrees with the hard equation of state without momentum dependence. 
However, the onset of squeeze-out in CASCADE mode is shifted towards larger impact parameters 
in comparison with the hard equation of state. Most importantly, the inclusion of the momentum 
dependence results in a drastic increase of the squeeze-out ratio. The lines are inserted to guide 
the eye. 
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